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A dynamic network in a dynamic population: 
asymptotic properties 

Tom Britton* Mathias Lindholm^and Tatyana Turova* 

January 18, 2013 

Abstract 



We derive asymptotic properties for a stochastic dynamic network model in a 
stochastic dynamic population. In the model, nodes give birth to new nodes until 
they die, each node being equipped with a social index given at birth. During the 
life of a node it creates edges to other nodes, nodes with high social index at higher 
rate, and edges disappear randomly in time. For this model we derive criterion for 
when a giant connected component exists after the process has evolved for a long 
period of time, assuming the node population grows to infinity. We also obtain an 
explicit expression for the degree correlation p (of neighbouring nodes) which shows 
that p is always positive irrespective of parameter values in one of the two treated 
submodels, and may be either positive or negative in the other model, depending 
QO . on the parameters. 
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1 Introduction 



The models of dynamical graphs are denned by the rules of attachment and deletion of 
vertices and edges chosen to fit particular processes in nature. Intensive study in this area 
began in the end of the 20th century when a number of models were formulated primarily 
in the physics literature (see, e.g., Barabasi et al. (1999), Callaway et al. (2001), Malyshev 
(1998)). The important paper by Bollobas, Janson and Riordan (2007) provided a unified 
approach based on branching processes to many known models of random networks (see 
the reference list in Bollobas et al. (2007)) 

In the present paper we continue the study of a Markovian model describing a random 
time-dynamic network in a random time-dynamic population. The model was originally 
defined by Britton and Lindholm (2010) extending an earlier model of Turova (2003), 
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which in turn was derived from a general model of Malyshev (1998). The population 
process is a Markovian linear birth-and-death process in which individuals are assigned 
random i.i.d. social indices. Given the population process, a Markovian network is de- 
fined where edges between individuals appear and disappear in such a way that "social" 
individuals tend to have more neighbours. Notice also, that it was shown by Turova 
(2002) that the model without social index and also without deletion of nodes includes 
as a subcase yet another model studied by Callaway et al. (2001). 

We study asymptotic properties of the network: what properties will the network (and 
population) have after having evolved for a long time assuming the size of the population 
has grown large. One can check that a snap-shot of the limiting network of Britton and 
Lindholm (2010) falls into the general class of inhomogeneous random graphs introduced 
by Bollobas, Janson and Riordan (2007). In particular, one can consult Britton and 
Lindholm (2010) as well as Bollobas, Janson and Riordan (2007) for the age, the type 
and the degree distribution in this model. 

In the present paper we shall make essential use of the theory of Bollobas, Janson 
and Riordan (2007) to study the phase transition in the model by Britton and Lindholm 
(2010). More precisely, we will determine the critical values which separate the set of pa- 
rameters under which the model with a high probability has a giant connected component 
(i.e. of order of the entire network), and the area of parameters which do not produce 
a giant component. This extends earlier results of Turova (2007) for the model without 
social index. 

Besides global properties of the network, e.g. phase transition, we shall study here a 
characteristic, which tells more about the local structure, namely the degree correlation, 
or the mixing coefficient as it was introduced by Callaway et al. (2001). Inspiring nu- 
merical results on the mixing coefficient in different empirical networks were presented 
and analyzed already by Newman (2002). Recently Bollobas, Janson and Riordan (2011) 
derived formulae for the mixing coefficient in terms of small subgraphs counts for a rather 
general graph model, which includes, in particular, the inhomogeneous graphs. 

According to empirical results provided by Newman (2002) real life networks appear to 
be of two classes: assortative, when the mixing coefficient is positive, these are primarily 
the social networks, or disassortative, when the mixing coefficient is negative, which is a 
feature of technological or biological networks. Still, most most of the models tend to have 
positive degree correlation/mixing coefficient. However, Newman (2002) could construct 
an example of a network model where the mixing coefficient changes sign depending on 
the parameters of the model. Lately Bollobas, Janson and Riordan (2011) gave another 
recipe to construct a network with a negative mixing coefficient. Therefore it is a challenge 
to look for an example of a somewhat "naturally grown" random network which possesses 
the property of enabling both positive and negative mixing coefficient depending on the 
parameters. As it is pointed out by Newman (2002) the property of assortativity affects 
qualitatively the sharpness of the phase transition with respect to the size of the giant 
component: the phase transition is sharper for the dissasortative networks. 

Here we provide an example of a dynamic (growing) social network with a degree 
correlation (mixing coefficient) that may be either negative or positive depending on 
model parameters. We derive an explicit formula for the degree correlation, showing 
the dependence of the sign of the mixing coefficient on the parameters. Notice, that 
the formula for the mixing coefficient for the model of Callaway et al. (2001) derived by 
Newman (2002), is a subcase of the formula we prove here. 



2 



Although our model is a subcase of the general model studied in Bollobas, Janson and 
Riordan (2011), and hence the general formula from this paper should be applicable for 
our model as well, we find a somewhat more direct way to get a formula for the degree 
correlation in our case. We use the invariant measures for the random walk associated with 
the graph. It is worth noting that our method is not restricted to the particular model we 
study here, but would work as well for a class of inhomogeneous random graphs. However 
this will be a subject of a separate study. 

2 The Markovian random network in a Markovian 
dynamic population 

Below we define the Markovian random network in a Markovian dynamic population, 
originally defined in Britton and Lindholm (2010). There are two versions of the model, 
referred to as the Uniform (U) version and the Preferential (P) version. In what follows 
a network denotes a finite set of nodes (the population) together with undirected edges 
connecting pairs of nodes. Nodes that are directly connected by an edge are called neigh- 
bours (in sociological applications nodes correspond to individuals, edges to some type 
of friendship, and neighbours are referred to as friends). The model is dynamic in the 
double sense that nodes are born and may die, and the same applies to edges. 

2.1 The model and its two versions 

We first define the node population dynamics. Let Y(t) denote the number of nodes alive 
at time t, and assume that Y(0) = 1. While alive, each node gives birth to new nodes at 
the constant rate A and each node lives for an exponentially distributed time having mean 
l//i (denoted Exp (//)), so each node dies at the rate \x. We assume that A > fi implying 
that the expected number of births during a life-span is larger than 1. This means that the 
node population is modelled by a Markovian super critical branching process. Additional 
to this, each node i is at birth given a random "social index" Si having distribution F s 
on M + , independent and identically distributed for different nodes. We will throughout 
assume that S has finite mean EfS] := fi$ < oo. 

There are two versions of the model for births and deaths of edges, both being Marko- 
vian given the node population and their social indices. In both versions, nodes are 
isolated at birth, i.e. have no neighbours. During the life of a node i, having social index 
Si = s, it creates new neighbour edges at rate as, also the same for both models. The dif- 
ference between the two versions lie in how neighbour nodes are selected. In the uniform 
(U) version the neighbour node is chosen uniformly at random among all living nodes. In 
the preferential (P) version of the model the neighbour node is instead chosen at random 
among all living nodes with probabilities proportional to their social index. Finally, in 
both versions each edge is removed, independently of everything else, at the rate f3. If a 
node dies all edges connected to the node in question are removed. We emphasize that a 
node gets new neighbours in two ways: it "creates" new neighbour edges itself but may 
also be selected as neighbour of another nodes that has created a neighbour edge. 
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2.2 Comments on the model 

The model has four parameters: the birth and death rates of nodes, A and \x respectively, 
and the death rate /3 of edges and a which is related to the birth rate of new edges. 
Beside these four parameters there is the distribution F$ for the social indices {Si} (in 
what follows we let f s denote the corresponding density function). 

The P- version of the model is inspired by the preferential attachment model [Ij, still 
being different in that here the probability of receiving edges is determined at birth 
whereas it is determined by random events during life in the preferential attachment 
model. In both models there is also an age-factor in the sense that older nodes tend to 
have more edges. 

Some submodels are worth mentioning. One is where there is no node-heterogeneity 
and nodes have the same social index S = s (for example set to 1 without loss of gener- 
ality). Loosely speaking, allowing node-heterogeneity makes it possibly to achieve degree 
distributions having heavy tails: if Fg is heavy tailed there will be some nodes with very 
high social indices that hence have a large number of neighbours (cf. the next section). 
The case where /i = and S = 1 has been studied by Turova (2002) and (2007) who 
derive more results for this submodel. 

2.3 Known results of the model 

The model allows for multiple edges and self-loops which usually does not make sense in 
applications. However, it was shown by Britton and Lindholm (2010) that the proportion 
of such edges is asymptotically negligible as long as E[S] < oo for the U- version, and as 
long as E[S 2 } < oo for the P- version of the model. As a consequence the network will 
then have identical properties if loops and multiple edges are ignored or not allowed. 

In the rest of the paper we will only consider the case where the node population 
tends to infinity, i.e. we condition on the event B : = {u; Y(t) — > oo}, which has positive 
probability since it was assumed that A > fi. The first two results follow from standard 
branching process theory: 

Asymptotic population size: As t — > oo 



Stable age distribution: As t — > oo the age A of a randomly selected living node will 
satisfy A ~ Exp(A) having mean 1/A, see e.g. Section 3.4 in Haccou et al. (2005). 

Social index distribution: Since the social indices are defined to be i.i.d. having distri- 
bution Fs, and they only affect the occurrence of edges, the social index of a randomly 
selected living node will have the same social index distribution F$. 

Asymptotic degree distribution: As t — > oo the degree (i.e., the number of neighbours) 
D of a randomly selected living individual has a mixed Poisson distribution in both 
versions of the model (Britton and Lindholm, 2010). For the U- version it is 



where A ~ Exp(A) and S ~ Fs are independent. This simply means that, conditional 
on A = a and S = s, the degree is Poisson distributed with mean parameter equal to 



Y(t) ~ e 



,t(A-/x) 




(1) 
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a(s + fis) (l — e (^+^) a ) + /i). The mean and variance of this distribution are given 
by 

4Aa 2 2 2a 2 

+ (A + (3 + /i) 2 (A + 2(/3 + n)f s + (A + (3 + li)(\ + 2(/3 + //)) [ J ' 



For the P-version of the model the degree distribution satisfies 

(2) 



, p , ( 2aS (1 - e-W+ri*) 
D {p) ~ MixPo I i '- 



where A ~ Exp(A) and S ~ Fs are independent. The mean and variance are given by 

4Aa 2 2 8a 2 

+ (A + (3 + ^) 2 (A + 2(/3 + ii)f s + (A + (3 + /i)(A + 2(/3 + fx)) [ J ' 

i.e. the same mean degree but a larger variance. For both versions it is seen that the 
variance of the degree distribution increases with the variance of the social index distri- 
bution, so having a heavy-tailed social index distribution Fs implies that also the degree 
distribution will be heavy-tailed. 



3 Results 

In what follows we derive some further properties of the model, still assuming t — > oo 
and conditioning on that the node population grows beyond all limits. We do this for 
the U-version and only present the results for the P-version which is obtained in a similar 
way. 



3.1 Type-distribution of neighbour nodes 

Assume t to be large and consider a randomly selected individual alive at this time 
having age a and social index s. We derive the distribution of the number of (a', s')- 
neighbours this (a, s)-node has. The number of neighbours of any type (a', s') will be 
Poisson distributed, so we first compute the mean of the distribution and then divide by 
the total mean thus obtaining the type-distribution. 

How many neighbours of type (a', s') does the (a, s)-individual have? Any such edge 
must have been created at some time r e (t — a A a', t) ('A' denotes minimum: both 
individuals must have been born). At such r our (a, s)-individual creates such an edge at 
rate asAe~ A ( T ~^~ a '^/ s (s / ), because it creates edges at rate as and the second part denotes 
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the fraction of nodes having age r — (t — a') (at r) and type s'. Similarly, it receives edges 
at the rate as'\e^ x( - T ^^ a ^ fs(s') because there are F(r)Ae _A ^ r_ ^ _a "fs{s') having the 
"correct" age and type, each of them creates edges at rate as' and with probability 
\/Y{t) it reaches our (a, s)-individual. Any such edge created at r remains at t if and 
only if both the (a', s')-individual and the edge survive until t (we ave already conditioned 
on that our (a, s)-individual lives at t). The probability for this equals e~^ + ^^ T \ The 
total expected number of edges with (a', s')-indivi duals hence equals 

{a', s'\a, s)da'ds' = da'ds 1 f a{s + s') Ae-^'-^/s^'K^^W. (3) 

J t—aAa' 

p(\—/3—ii)aAa' i 

= a(s + s')f s {s')\e- Xa '— x da'ds'. (4) 

It is easy to check from this that the expected total number of neighbours of any type of 
the (a, s)-node, the integral of this quantity over all a' and s', equals 

(u) _ a{s + fls) h _ -03+ M )a 



m y M-,-a,s 



(1 - e -^ a ) . (5) 



This also confirms that the unconditional degree distribution has the mixed Poisson dis- 
tribution specified in Equation (JTJ). 

As a consequence, for the U-version of the model, the type distribution of the neigh- 
bours of an individual of type (a, s), which equals the expected number of (a 1 , s')-types 
divided by the total expected number of types, becomes 

f (U \a',s'\a,s) 



mW(a',s'\ 


a, s) 


m( f/ )(-,-| 


a, s) 



s + fis (A - $ - fi) (1 - e-(P+ri a ) 

We hence have conditional independence: f^ u \a', s'\a, s) = f^ u \a'\a)f^ u \s'\s), where 

(B + u)\e~ Xa ' (fi^P-v)^' - i) 
1 [ala) (X-p-Ml-e-V+i*) ' [) 

fV\s'\s) = + (7) 

Admittedly, the notation here, and in what follows, is a bit sloppy in that f <yU \a'\a) 
and p u '{s'\s) really denote different density functions, as indicated by the difference in 
arguments. 

Using similar arguments, but for the P- version of the model, we get f( p \a', s'\a, s) = 
f( p \a'\a)f( p \s'\s) where 

f( p )(a'\a) = [ ' N [ - -J- (8) 

J 1 1 ' (A-/3- /U )(l- e -(^) ' [) 

/ W Ma) = £W) > (9) 
i.e. the same density with respect to age but different with respect to social index. 
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3.2 Phase transition 



Here we shall apply the results of Bollobas, Janson and Riordan (2007) to the introduced 
above models in the limit as t — > oo. Therefore first we shall put our model into a general 
setup of the theory (and notation) of inhomogeneous random graphs of Bollobas, Janson 
and Riordan (2007). 

Let us enumerate the individuals at time t by % — 1, . . . , n :— Y(t). Let Xi = (Ai, Si) 
denote, correspondingly, the age and social index of individual i. Notice, that this means 
that this individual was born at time t — Ai (and remains alive at time t). The random 
variables Xi are independent for different i. For a randomly selected i, as we argued 
above, the distribution of the social index follows the original social index distribution 
and Ai is Exp(7); denote the corresponding measures by /xi and p, 2 . 



Hi(da) = \e Xa da, 
/i 2 (ds) = f s (s)ds. 



(10) 



We shall call Xi the type of individual i. Given the types Xi = (ai, Si) and Xj = (a,j, Sj) of 
individuals i and j, we derive the probability Pij(n) that they are connected. From above 
we know that the expected number of (a', s')-indivi duals our i is connected to equals 
m(a' , s'\ai, Si)da'ds' defined by 01]) for the U-version (and similarly for the P-version). 
In total, if Y(t) = n, the expected total number of individuals of type (a', s') equals 
n\e- Xa 'da'f s (s')ds'. This reasoning implies that the probability that i is connected to j 
(of type (a,, Sj)) equals the ratio of these two expressions: 

Pij(n) = -Ki(ai,aj) K 2 (si,Sj), (11) 



where 



and 



e 



(X-/3-fi)(a i Aa j ) _ j 



Kxicii, a 



k 2 (s^ Sj) = a(s + s') in U-version, 



K 2 {Si, s 



K 2 { s ii s j) = 2ass' in P-version. 



Let fj, = x /i 2 , which by ffTUj) is a probability measure on S = R^, and let 

K = K,(Xi,Xj) = Ki(a i) aj)K 2 (s i) Sj). 

Given the sequence xi, . . . , x n , we let G v (n, k) be the random graph on {1, . . . , n}, such 
that any two vertices % and j are connected by an edge, independently of the others 
vertices, with probability given by (|lip . which we write as 

Pij{n) = min{fi;(xi, Xj)/n, 1}. 

Hence, conditionally on Y(t) = n the graph G v (n, k) describes our model. On the other 
hand, G v (n, k) with k as above satisfies the definition of inhomogeneous random graph 
from Bollobas, Janson and Riordan (2007), and thus we can apply some known results 
here. 



7 



Recall first the fundamental result on Phase Transitions in the inhomogeneous random 
graphs. Define 

T K f(x)= / K{x,y)f{y)d/i{y), 



and 

||Tj=sup{||T K /|| 2 :/>0,||/|| 2 <l}. 

Then, by Theorem 3.1 from [3], the largest connected component of G v (n, k), denoted 
Ci(G v (n, k)), satisfies 



-C 1 (G v (n, k)) 4 p K = [ p{x)dp( 
n Js 



s x), as n — > oo, 

Js 

where p K > if and only if \\T K \\ > 1. 

This gives as well for our model that conditionally on the number Y(t) of individuals in 
the population at time t, we have with a high probability approximately p K Y(t) individuals 
connected; and thus this holds unconditionally as well. The function p K (x) is described 
more precisely in the following theorem. 

Theorem 3.1 ([3]. Theorem 6.1) Suppose that k is the kernel on (S,p), that k G L x , 
and 



K(x,y)dfi(y) < oo 
for every x 6 S. Then p K (x) is the maximal solution to 

f(x) = l-exp{-TJ(x)}. (12) 

Furthermore: 

(i) If \\T K \\ < 1 then p K {x) = for every x, and KW\) has only the zero solution. 

(ii) If 1 < \\T K \\ < oo then p K (x) > on a set of a positive measure. If, in addition, k 
is irreducible, i.e., if 

ACS and k = a.e. on A x (S\A) implies p(A) or p(S \ A) — 0, 

then p K (x) > for a.e. x, and p K (x) is the only non-zero solution of (Q] 



First we shall derive conditions for p K (x) > 0. By the cited results this amounts to 
computing \\T K \\. Unfortunately, there is no general formula to simply compute the norm 
||T K ||. Therefore we shall use a criteria for the condition when ||T K || > 1, which is what 
we need here. 

Proposition 3.2 /J3J/, Proposition 17.2] For k > 1 let 

a(k) := - I k(x ,xi)k(xi,x 2 ) . . . K(x k -i,x k )dp,(x Q ) . . .dp(x k ). 

I Jsk+i 

Then \\T K \\ > 1 if and only if a(k) — > oo as k — )■ oo. 
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It follows from definition ffTUl) that for our model 



where for i = 1, 2, 



OLi{k) 



R 



fc+i 



a(fc) = -oti(k)a 2 {k), (13) 



«i(wo, W2) • • • Ki{uk~i, u k )dni(u ) . . . d/ii(u k ). (14) 



First we shall compute a2{k). We start with the U- version in the next Proposition. 
Proposition 3.3 Introduce the matrix 

E[S] 1 



A 



E[S 2 } E[S] 



Then for all k > 3 



a u 2 {k) =a k \A k - 1 



2E[S] 
[E[S}) 2 + E[S 2 ] 



(15) 
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where for a matrix B we denote {B}- the corresponding entry 



Proof. Let So, Si, . . . , be i.i.d. copies of the random variable S. Then by the definition 



a$(k) = a K E 



\[{Si-i + Si 



t=i 



=: a h f(k), 



(16) 



where 



Denote also 



f(k) = E 



HiSi-i + S, 



g{k) = E 



' k \ 
l[(S-i + ^) S* 

\i=l J 



Then we recursively derive 



and 

for k > 1, with 
Hence, for all > 1 



/(*) = 0(*-l)+/(*-l)E[£], 

^(*) = ^-l)E[5]+/(*-l)E[5 a ], 
/(1)=2E[5], s(l) = (E[S]) 2 + E[S 2 ] 



(17) 



/(*) 



.4 



^-1) 

and formula (1151) follows from here, (fPTj) and ( 1T61) 



1 ! 9(1) 



□ 
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Corollary 3.4 One has 



lim (a^(k)) 1/k = a ( E[S] + y/BS^] 



(18) 



In the case of the P-version, when «2 = ac% h is straightforward to derive from defini- 
tion flUD that 



a. 



\k) = (2a) k (E[S}Y(E[S 2 ]) 



which yields 



Define now 



lim {a%{k)) 1/k = 2aE[S 2 



k— >oo 



c cr (A,/i + 0) := sup{x > :^2x k a x {k) < oo}. 



(19) 
(20) 

(21) 



k=2 



Then due to the definitions (ITS]) and (fill) together with asymptotics (ITS]) and (1211 we 
have the following criteria which holds for both cases, P-version and U-version. 



Corollary 3.5 // 

then a(k) — > oo; if 



lim (a 2 (k)) 1/k > c cr (A,/i + /3) 
lim (a 2 (£;)) 1/fe <c cr (A,/i + /3) 



k— >oo 



i/ien a(/c) < oo. 

The value c cr (A,/i + /3) is known from [11] . Let us record this result here. Write first 



oii(k) 



R 



fc+i 



Ki(ao, ai) • . . . • Ki(afc_i, afc)Ae Aa °<iao . . . Ae Xak dak 



Ki 



Xq Xi 

T'T 



-<*fc-l Xk 



where X , . . . , are z.i.d. random variables with a common Exp (l)-distribution. Then 
by the results of [12] (see [12] formula (1.7), or consult [3]) the value c cr (X,fi + (3) is the 
smallest positive root of 



oo / \ n 1 n 1 



(22) 



Finally, combining Corollary 13.51 and Proposition 13.21 together with asymptotics (ITBl 
and (1201) we can derive, apart from the critical case, the necessary and sufficient conditions 
for the existence of the giant component. 



Corollary 3.6 Denote 
and let 



u _ u p _ P 



a 



E[S] + ^EW] 

c cr (X,fi + (3) 



and R p (a, /3, A, 7, S) :- 



2aE[S 2 



(23) 
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Then there exists a giant component in the U -model if R u (a,(3, A, 7, S) > 1, and there is 
no giant component if R u (a, (3, A, 7, S) < 1. 

Similarly, there exists a giant component in the P -model if R p (a, (3, A, 7, S) > 1, and 
there is no giant component if R p (a, (3, A, 7, S) < 1. 

Note also that although we define the function c cr (X,u) for u > 0, the model with 
p + (3 = is treated in a very similar way, in fact when also S = 1, this model becomes 
the one studied by Callaway et al. (2001). Then one defines c cr (A,0) as the critical 
parameter, above which there is a giant component. It is known (see Callaway et al. 
(2001)) that c cr (A,0) = 1/4, and moreover function c cr (A,w) is continuous at u = (see 
Bollobas et al. (2007) and Turova (2007a)). 

In general, there is no closed form for the function c cr , but one can mention some 
qualitative properties of this function. First of all, simply from the model it follows that 
c cr (A,/i + (3) is decreasing in A and increasing in p + fl, i.e., both in p and (3. For any 
u > denote c{u) the smallest positive root of the rescaled function (|22|) 

00 1 

H u (x) = i + y2(-irx n - 

n=l 

Then 

c cr (A,/i + /3) = (jm + P)c 



n 

1=1 



1 + (l- l)u 




3.3 Degree correlation 

We now derive an expression for the degree correlation p (also known as "mixing coeffi- 
cient") of the network when it has grown and reached its stationary phase. That is, after 
a long time t we take a snap-shot of the network and compute the degree correlation p, 
which is defined as the correlation of the degrees of the two adjacent nodes of a randomly 
selected edge in the network. We derive p in two different ways. In the first method we 
first pick a random node in the snap-shot network and then perform a random walk on 
this fixed network. When this random walk has reached stationarity we consider the node 
in one step as the "first" node and then pick the "second" node randomly among the 
neighbouring nodes, and let the edge between these two nodes be our randomly selected 
edge. In the second method we pick our "first" node randomly in a size-biased way in the 
snap-shot network (such that the node is adjacent to a randomly selected edge) and the 
"second" node randomly among the neighbours of this individual (the size-biased distri- 
bution of a non-negative random variable X with density/pdf fx{x) and finite mean px 
has density/pdf xfx{x) / px - below we denote such a random variable by X). As will be 
seen, the two methods of picking neighbouring nodes give the same result. The degree 
correlation is then obtained by computing the degree correlation of the two neighbouring 
nodes. 



3.3.1 The type-distribution obtained using random walk method 

We start with the first method where we first pick a node at random among all living nodes. 
From the results of Section 12.31 this means the node has exponential age distribution and 
independent social index distribution given by the original distribution. More precisely 
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the type distribution of a randomly selected node equals Xe~ Xa dafs(s)ds. The neighbours 
of this node, conditional on its type, all have distribution f(a', s'\a, s), defined in Section 
3.1, where its form depends on whether we consider the U- version or the P- version of 
the model. Iterating this procedure gives the type-distribution of nodes after additional 
steps in the random walk. Eventually this random walk reaches stationarity and the 
distribution foo(a', s') of the type of node is then given by the solution to the functional 
equation 

foo(a',s')= / / f 00 (a,s)f(a',s'\s,a)dsda. 



And, since f(a',s'\s,a) = f(a'\a)f(s'\s) in both versions of the model, it follows that 
the age and social index distributions are independent in the stationary distribution as 
well: foo^a'.s') = foc{a')foo{s')- The corresponding densities are hence the solutions to 

/ooO') = J foo{a)f{a'\a)da, (24) 

foods') = J fMfa'Wds. (25) 

The functions f(a'\a) and f(s'\s) were defined in Section 3.1, f(a'\a) being the same 
for the two version and f(s'\s) being different for the U- version and the P- version. The 
(unique) solution /oo(a') to ([21]) is given by 



/oo(a') = A (l + -^-\ e~" (l - e-V+M) . 



(26) 



As for the solution to ( 1251) . its solution depends on which version of the model we consider, 
i.e. whether we use f^ p \s'\s) or f( u '(s'\s). The (unique) solutions f£\s') and /£V) 
are given by 

f£\s') = s>f s (s>)/^ (27) 

and 

/<£°0O = "-^fsis'). (28) 

The P-version hence has the size-biased version of the social index distribution whereas 
the U-version has a mixture of the original and the size-biased social index distribution. 

The stationary distribution defined in (|26|) and (127|) for the P-version, and (1261) and (1281) 
for the U-version, is hence the distribution of the "first" node of a randomly selected edge. 
Given the type (a', s') of this individual, a randomly selected neighbour of this individual 
has the previously defined type-distribution f(a,s\a',s') = f(a\a')f(s\s r ) (where f(s\s'), 
but not f(a\a') depends on which model version we consider). 



3.3.2 The type-distribution obtained using the size-biased method 

We now derive the type distribution of the nodes of a randomly selected edge using the 
second method described above, which will be seen to give the same solution. We do this 
for the P-version of the model. 

We know from before (see Equation |2]) that the degree D of a randomly selected node 
in the snap-shot network is mixed Poisson with (random) parameter 

2a5(l-e- (/3+ ^ )A )/(/3 + /i), (29) 
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where S has the original social index distribution fg and A is independent Exp(A). We 
now instead pick our "first" node as a node of a randomly selected edge. As a consequence, 
this node will have the size-biased version D of this random variable (because a node of 
degree k has probability proportional to k of being selected). The size-biased distribution 
of mixed Poisson is also mixed Poisson, but with the size-biased version of the random 
parameter. And, the size-biased distribution of the random parameter (J2HJ) equals 




By this tilde- notation we mean that the type (A*, S*) of the "first" node will have social 
index distribution S (the size biased distribution of the original density, having density 
sfs(s)/fjLs) which hence agrees with (J27|) . The age distribution A* will be independent 

and have the distribution defined by 1 — e~^ + ^ A * being distributed like ^1 — e~(^+^) A j, 

where A ~ Exp(A). But, since A ~ Exp(A) it follows that 

i _ e -v+»)A „ Beta (Y A 



+ v 

' From this it follows that the size biased distribution 



(l— ^)~Beta(2, 



So 1 - e~ (/3+M)A * has this distribution, implying that e" (/3+M)A * ~ Beta (A/(/3 + fi), 2). 
From this it follows after some straightforward calculations that A* has exactly the density 
given in Equation (l2T)j) . As in the previous sub-section, given the type (a', s') of the "first" 
node, a randomly selected neighbour of this node has the type-distribution f(a, s\af, s') = 
f{a\a')f{s\s'). 

3.3.3 The degree correlation p 

In the two previous sub-sections we have derived the type-distribution of the two nodes of 
a randomly selected edge: The type (A*, S*) of the "first" node has density /oo(a')/°o( s ') 
where foo(a') was defined in ( )26l) for both versions of the model, and foo(s') was given 
by ( 12 7p for the P- version (and not explicit for the U- version). Given the type (a', s') of 
the first node, the second node has type-distribution f(a,s\a',s r ) = f \a\a') f (s\s') , where 
f(a\a') and f(s\s') were defined in - fl9]). We now derive the degree correlation between 
two such randomly selected nodes. 

In order to compute the degree correlation of the nodes adjacent to the randomly 
selected edge we condition on the types. Given the types we know from before that 
the degree distribution is Poisson: an individual of type (a, s) has a Poisson number of 
neighbours and the mean g(a, s) = gi(a)g2(s) equals 

l- e -^+^ a \ gf\ s ) =a(s + ns) for the U-version, 

p + fM I <?2 ( s ) = 2«s for the P- version. 

(30) 

We now compute the degree correlation using these results. Let Di denote the degree 
of the "first" node and D 2 to its selected neighbour. In order to compute the degree 
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correlation p(D 1 ,D 2 ) we need to compute the following expectations: E[Di] = E[£) 2 ], 
E[Df\ = E[£)|] and E[DiD2] (the first as well as the second moment of D\ and D 2 are the 
same since we have reached stationarity). We do this excluding the edge connecting the 
two nodes (subtracting the value 1 everywhere does not change the degree correlation). 
The conditional degree distributions remain Poisson with the same means because if 
X ~ Po(/3), then the conditional distribution of X — 1, conditional on that X > 1, is also 
Po(/3). We hence get 



E[Di] = E[D 2 ) = J gx{a)f oa {a)da J g2(s)f 00 (s)ds, 
E[D\] = E[Dl\ = J J{g 1 {a)g 2 {s) + gl{a)gl{s))f 00 {a)f 00 {s)dads ) 
E[D 1 D 2 } = J j g 1 (a)g l (a')f 00 (a)f(a'\a)dada' J Jg2(s)g 2 (s')f O0 (s)f(s'\s)dsds'. 
Given these expressions, the degree correlation is given by 

p{Du ° 2) ~ E{Dl) - (E(D 1 )) 2 ~ V[A] • (31) 

The perhaps most important questions is to learn if p, or equivalently the covariance 
C, is positive or negative. For this it is sufficient to compute the numerator of f pJTj) . 
Using /oo(a, s) (defined by ([25]) . (127]). and f[2"51) ) and (15U1) in the expressions above for the 
P-version and U-version, separately, standard but tedious calculations reveal that 



E[-D} ] = E[Z^ 



A + 2(/3 + p) E[5] 



(p) n (Ph_ 5A + 6Q0 + //) /2«E[5 2 



(A + + n)(\ + 2(/3 + /i))(A + 3(/3 + fi)) V E[5] 
as well as 

EfD^l - EfD^l - 2 Q (E[g 2 ]+3(E[S]) 2 ) 

This gives us for the P-version 

C ' P,(g " = (A + g + „)(A + 2(/I „ OT A + 3(g + gj < 2aE(S2) ) 2 • (32 > 

Since all parameters are positive we conclude that for the P-version of the model the 
covariance, and hence also the degree correlation p, is always positive. This is true 
irrespective of the model parameters and choices of social index distribution. 
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For the U- version the picture is different. Introduce 7 = /3 + p and compute 



-a 2 



1 fE[S 2 ] + 3(E[S]) 2 ^ 2 



(A + 2 7 ) 2 V E[5] 



« 2 r(E ^ )2 (67 + 5A) (A + 2 7 ) 2 / Egl + ^ _ f JgL + 3 



(A + 2 7 ) 2V ' \ (A + 7 )(A + 3 7 ) V(E[5]) 2 ; V(E[5]) 



(34) 



Let us also denote 

A 2 



(A + 7 )(A + 3 7 )' 

Then we can rewrite fl33l) as follows 

C {U) (D U D 2 ) = j^—jiES) 2 ( 2(a + 4) f Jfl + 1 )- ( Jfl + 3 V j (3-V) 



(a + 2 7 ) 2 ^ ' v v y V(E[5]) 2 y V(E[.V|) 

(E^-fS) 2 + ^ + D-d- 2a) 



(A + 2 7 ) 2V ' V V(E[5]) 2 ; (E[S)) 
This together with the facts that (fjfjp- > 1 and < a < 1, yields that if 

— f— ^ < 1 + a + Va 2 + 4a 

then C^(Di,D 2 ) > 0, and hence also that pM > 0. If 

E[S 21 



;e[5]) s 



+ Va 2 + 4a 



then = 0, and otherwise, p^ u ' < 0. 

We conclude, that for a fixed value of E[S] and the other model parameters, increasing 
the variance of S from (constant S) to very large variance, allows us to pass through the 
assortative regime, to neutral (no assortativity) and then to disassortative. One can also 
observe that increasing the variance (i.e. EfS* 2 ]) also increase R u and R p respectively (see 
(1231) ) which thus might cross the critical value of 1 leading to percolation. This second 
effect, that disassortativity facilitates percolation, was also reported by Newman (2002) 
for a different model. 



4 Discussion 

In the present paper we studied properties of a stochastic dynamic network in which 
nodes are born and die randomly in time, and while alive pair of nodes are connected 
and disconnected by edges randomly in time. The rate of connecting to other nodes 
depended on the social index, given at birth, of the node. For this model we derived 
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limiting properties of a snapshot of the network after a long time, assuming the node 
population grew large. 

The three main results were: the type distribution of neighbour nodes (Section 3.1), a 
criterion for when phase transition occurs, above which the network has a giant component 
(Section 3.2) and an expression for the degree correlation of connected edges (Section 3.3). 
One of the most interesting features of the model is that U-version of the model (in which 
neighbouring nodes are selected uniformly among all living nodes) may have either positive 
or negative degree correlation depending on the numerical values of the model parameters 
(birth and death rates of nodes and edge) and the social index distribution. 

Together with previous results of the model (Britton and Lindholm, (2010)) the most 
important limiting local properties of the network are hence known by now, as well as 
whether or not the network has a giant component. Other global properties such as the 
diameter of the network remains to be analysed. Another interesting class of problems to 
study for the present model would be to study dynamic properties of the network. One 
such problem would be to study limiting properties of some process taking place "on" 
the network. For example, can an epidemic, forest fire or similar persist on the network 
forever or will it die out? What happens if a lightning process kills randomly selected 
nodes and all their neighbours at a constant rate. 

The model can of course also be generalized in several ways to make it more applicable 
for certain situations. Nodes could be of different types with different attachment rates 
depending on the types of the nodes in question, for example male or female if mimicking 
a sexual network. Similarly, edges could be of different types, perhaps reflecting the 
"degree" of acquaintance which for example could affect the transmission probability of 
an epidemic taking place on the network. Still, the current model does capture some 
important properties of a real world network in the nodes as well as edges are created 
and cease to exist. Also, depending on the model parameters, the model is quite flexible 
in producing different degree distributions and degree correlations (the network has no 
clustering asymptotically). 
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